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Abstract 

In this paper we statistically characterise the bio-aerosol background in an urban environment. To do this 
we measure concentration levels of naturally occurring microbiological material in the atmosphere over a 
two month period. Naturally occurring bioaerosols can be considered as noise, as they mask the presence 
of signals coming from biological material of interest (such as an intentionally released biological agent). 
Analysis of this 'biobackground' was undertaken in the 1-10/xm size range and a 3-9% contribution was 
found to be biological in origin - values which are in good agreement with other studies reported in the 
literature. A model based on the physics of turbulent mixing and dispersion was developed and validated 
against this analysis. The Gamma distribution (the basis of our model) is shown to comply with the scaling 
laws of the concentration moments of our data, which enables us to universally characterise both biological 
and non-biological material in the atmosphere. An application of this model is proposed to build a framework 
for the development of novel algorithms for bio-aerosol detection and rapid characterisation. 
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1. Introduction 



scope of the current paper and the readers are 



It is well recognised that biogenic aerosols 
have a significant effect on human health, the 
environment , clim a te, a t mosphere and othe r 
phenomena (IPoschll (l2005l ): Ijaenicke et all (|2007|) : 
Ho and DuncanI (|2005f ): ISuffman c t all (|2010| ) 
and references therein). Numerous studies have 
been conducted that focus on the characterisation 
of naturally occurring microbiological material 
in the atmosph ere (referred in t his paper as 
biobackground) ( Bauer et al. ( 20081 ) V This large 
body of literature indicates a rich and diversified 
dependence of biobackground characteristics af- 
fected by geographical location, temporal, seasonal 
and daily variations as well as t he dependency 
on ap plied measuring techniques (|Ho and Spence 
( 19981 )). A comprehensive review of publications 
concerning bioaerosol studies is outside of the 



(Huffman et al. 


(12010 


iHarrison et al. 


J2004f 



Burrows et all ( 2009 ) 



Jaenicke et al.l 



20071) 



lAriva et al.l ()2009D and references therein). 

Biobackground is the primary limiting factor 
of bio-detection technologies used in operational 
applications (bio- security, defence syste i ns and 



counterterrorism) (jHo and DuncanI (|2005l ): iJensen 
(|2007l) ). 

Naturally occurring bioaerosols can be considered 
as noise, as they mask the presence of signals com- 
ing from biological material of interest (e.g. an in- 
tentionally released biological agent), hindering a 
biodetectors' capability (by affecting its sensitivity, 
selectivity, time of response or rate of false alarms) 
to detect airborne biological threatfl 

Selection of a specific technique for bioaerosol 
detection is a challenging task, whose outcome is 
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determined by both the intended ap phcation and 



the sc enario in which it is deployed (|Buteau et al 
( 2010( 1). It is unhkely that a single biodetection 



technique capable of detecting all biological agents 
will be developed in the foreseeable future. This 
necessitates the use of multiple, often complemen- 
tary technologies and data processing algorithms 
for reliable bio-threat detection. The consistent in- 
tegration of these technologies into a single opera- 
tional system is influenced by the critical issue of 
'data fusion' of multiple information sources, i.e. 
the rigorous comparison of the relative informa- 
tion gain that can be provided by a particular bio- 
detection technique. Apparently this issue cannot 
be thoroughly investigated and eliminated without 
detailed knowledge of the biobackground the pro- 
posed system is operational in. From this perspec- 
tive; understanding the characteristics, dependen- 
cies and parameters affecting the biobackground 
is a critical step in the development of enhanced 
biodetection capability operating over a broad spec- 
trum of applications and scenarios. 

Particles in the ambient environment (atmo- 
spheric surface layer) are in continuous chaotic mo- 
tion being advected by turbulent flow. Since most 
of these particles are very small (less than lOfim) 
their inverse effect on flow is negligible and their 
kinematics can be described by the models of 'pas- 
sive scalars' (i.e. a system of passive markers associ- 
ated with the i ndividual particles rand omly moving 
by turbulence) ( Falkovich et al. ( 200l[) ). Under this 
condition it is reasonable to assume that the funda- 
mental physical properties of turbulent mixing and 
not the physical properties of the particles (e.g. size, 
density etc.) would strongly emerge in the statis- 
tical characteristics of the biobackground. These 
characteristics are determined by the kinematics of 
the turbulent flow and are therefore independent of 
the origin or nature of the particles (i.e. biological 
or non-biological). 

An illustrative example of this approach 
has been recently re ported in our publication, 
Skvortsov et al. ( 2010f) . where we established scal- 
ing laws of time series for particle concentra- 
tion and found they com ply with specif ic mod - 
els of turbulent mixing (|Yee and ChanI 
Villermaux and DuplatI ( 2003 )). These results pro- 
vide the motivation and foundation for the current 
study, i.e. the application of the theoretical meth- 
ods of passive tracer dispersion in the atmospheric 
surface layer to deduce the universal statistics (scal- 
ing laws) of particles regardless of their origin. 



The current study focuses on a bio-detection 
framework which enables a consistent predic- 
tion of some of the statistical properties of par- 
ticle concentrations in the biobackground (mo- 
ments, probability distribution functions (PDFs), 
intermittency corrections, etc.) by employ- 
ing well-known propert i es of turbulent mix- 
• '" !e and Cha n (lOOyV^ ^ 
( 20031). iLebedev and Turitsvn (2004 ). ,Yee (.200E 



i ng (|Yee and Cha n* (TooT), Villermaux and DuplatI 

. n t " 

Skvortsov and Yed ( 20111 )'). The statistical frame- 
work presented was validated using a large data 
set of biological and non-biological particle concen- 
trations measured by an optical biodetector util- 
ising laser induced auto-fl uorescence cha r acteris - 
tics for living organ is ms (iHairston et ahl ( 1997 ) ; 
Ho and SpeiQ (|l998l) ; lBrosseau et all (|2000l )). We 
also discuss the applicability of this statistical ap- 
proach for the development of detection algorithms 
in operational systems. 

The framework was developed in two steps: first, 
we collected a large data set of physical charac- 
teristics of the ambient biobackground. The data 
were screened, validated, analysed and results com- 
pared to literature. This provided the quantifica- 
tion of a validated biobackground data set for an 
urban/industrial type of ambient environment. We 
then developed a mathematical model based on sta- 
tistical properties (such as moments) derived from 
analysis of time series of measured data, revealing 
analytical forms of concentration PDFs character- 
istic of the ambient biobackground. The derived 
generic PDFs (Probability Distribution Functions) 
are fundamental to the development of the next 
generation of biodetection algorithms. 

The paper is organized as follows: in the first part 
we outline the measuring methodology and present 
quantification results of the biobackground; in the 
second part we present the statistical model, its val- 
idation and the analytical PDF forms derived from 
our measured data. 



2. Experimental 

2.1. Instrumentation 

The atmospheric aerosol was measured and char- 
acterised by two instruments: Fluorescent Aerosol 
Particle System (FLAPS-III, Dycor Technologies) 
and Aerodynamic Particle Sizer (APS, TSI Model 
3321). Both have been used in a variety of aerosol 
research studies including the areas of environmen- 
tal and biobackground monitoring (jHuffman et al 
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( 20101 ) and references therein), hence only a brief 
description follows. 

FLAPS-III is an integrated optical biodetection 
system providing near real time (few seconds) 
point -detect i on of particles of biological ori- 
gin (iDvcoij (l2011al) ). It utilises laser induced 
auto-fluorescence intrinsic to li ving biological 



orgaii is ms ( Boulet et al. (Il996l): Hairston et al 



( 19971): Iffiil et all (Il999l): IHo and Spencd (| l99F 
Eng et all (Il989l): [Setlow and Setlowl (Il977 



Laflammc e t al. 120051) 7 lAgranovski and Ristovskil 
^005 ): Huffman et al.l (|2010l) ). The system is com- 
prised of two major components: FLAPS-III (TSI 
Model 3317) and particle concentrator (XMX/2A, 
Dycor Technologies). FLAPS-III measures side 
scattering (SS) and fluorescence (FL) light intensi- 
ties of a single particle. The fluorescence signal is 
measured in two wavelength bands corresponding 
to the emissions from NADH (FLl) and riboflavi n 
(FL2) ( "Hairsto n et al.l ()l997[ ): iBoulet et al.l (|l996[) : 
Hill et al. (199^) ). In conjunction with the XMX 
concentrator the system provides count rates of 
total (TAP) and fluorescent biological aerosol 
particles (FBAPs) in the l-lO/xm size range. The 
system does not allow identi fication or spe ciation 
of the detected bio-particles (jDvcoii (|2011bl )). De- 
tailed descriptions of the instrumen t and operating 



princi ple s are presente d elsew her e (iHairston et al 




IPinnick et al. ( 19981) : Agranovski et al 



FLAPS-III is currently considered as one of 
the benchmark biodetectors for use in defence 
applications. The system has been evaluated 
for accuracy, sensitivity and reliability in numer- 
ous laboratory and field studies using biological 
agents and simulants and has pro ved to be a re- 
liable a nd accura te bio-d et ector (IHo and Spence 



( 19981): lHoetal.1 (Il999l ): ISemler et all (l2004l): 



Huffman et al 



Agranovski and Ristovskil ( 2005[) ) 
(|2010l ) demonstrated the suitabilitv of FLAPS-III 
for field studies and biobackground characterisa- 
tion. 

APS is a time-of-flight optical spectrometer mea- 
suring number concentration and size characteris- 
tics of particles in the 0.5-20/ito size range. For 
the purposes of this study the APS measure- 
ment data was truncated to the l-lO/im size range 
to match data measured by FLAPS-III. Assum- 
ing particle sphericity and known density; sur- 
face, volume and mass characteristics of parti- 
cles can be estimated from the measured data. 
More detail about the APS instrument is pre- 



sented in (iBaron and Wilier (l2005l): ItoB (l201ll) : 
Peters and LeithI (|2003[) : iTsai et al.l (|2004l) ). 



2.2. Instruments settings and calibration 

APS and FLAPS-III were factory calibrated prior 
the commencement of the study and maintained 
during the measuring campaign according to their 
respective manuals. The operational settings for 
FLAPS-III recommended by the manufacturer were 
experimentally validated u sing aerosol ch a llenge s 



of known characteristi c s (iKanaani et aLl (|2(308l) ; 
Agranovski et al.] ( 2004 ): Huffman et al. r(|2010l)). 

The SS and FL signals measured by FLAPS-III 
are each gated into overall 32 bins according to 
the increasing SS and FL light intensities. Based 
on experimental results which were in agreement 
with the manufacturer's recommendation, bins [2- 
31] (where bin is the first bin) were selected for 
the measurement of fluorescent (FLl and FL2) sig- 
nals, and bins [0-31] for the measurement of SS 
signals. Settings for LED power and PMT gain 
were determined using non- fluorescent PSL spheres 
and NaCl aerosols following the method outlined 
in Brosseau et al 2000. FLAPS-III operated in a 
continuous sampling mode with the concentrator 
flowrate set at maximum (350-L/min) and aerosol 
stream flowrate IL/min. 

2.3. Fluorescent Biological Aerosol Particles 

Fluorescent biological aerosol particles (FBAPs) 
measured by FLAPS-III represent a subset of 
primary biological aerosol particles (PBAPs). 
PBAP can be defined as a particle that dis- 
cernibly is or was all or part of a liv i ng or - 
ganism dCabev et all (l2009f) : IPinnick et al!* (2009|); 

Ho and Spence (19981 ) 
FBAPs measured by 



.Huffman et al. 



(l20inl) ) 

demonstrated that the 
FLAPS-III technology are attributed mainly to 
the presence of biological material. The measured 
count of FBAPs may be considered as a conser- 
vative estimate of the PBAPs' a bundance in the 
ambient air (Huf fman et al. ( 2010l) ). 

Some particles, despite being of non-biological 



or igin, a lso exhibit fluorescence (|Agranovski et al 



(|2003a| ): ISivaorakasam et al.l (|2004r) ). however the 
contribution of these interferents to the overall FL 
counts measured by FLAPS-III for most ambient 
conditions and the size range of i nteres t (LTO^mx) 
can b e neglected (jPinnick et al. (l2009l) : iHo et al 
( 19991 ) and references therein). 



Since FLAPS-III does not provide direct quantifi- 
cation of aerosol concentration - only an integrated 
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count of particles detected in the sampled volume 
is available - it was used in parallel with the APS. 
The concentration of FBAP was then estimated 
from the total particle concentration measured by 
APS and the fraction of fluorescent particles in the 
total count was determined from FLAPS-III data 
(FL/SS ratio). 

3. Methods 

3.1. Air Sampling 

The measurements were conducted on the top of 
a four story building at 12m height with the in- 
strumentation housed inside an air-conditioned en- 
vironmental enclosure. The ambient air was sam- 
pled through vertical conductive sampling lines pro- 
truding about 1.5m outside of the enclosure's ceil- 
ing. Deposition losses of particles in sampling lines 
were assessed and found to be negligible. Sampling 
was done on a semi-continuous basis with the inter- 
ruptions due to instrument maintenance and avail- 
ability. Both APS and FLAPS-III operated side by 
side and provided measurements every 5 seconds. 

3.2. Measurement location 

The sampling site was located near Port Mel- 
bourne (-37°49'27.45", -Hl44°54'43.02"), approx- 
imately 0.5km South from the Yarra River; 5km 
West from the Melbourne CBD and 3km North 
of Port Philip Bay. The site is representative of 
an urban/industrial environment influenced by ma- 
rine sources during SE to SW wind conditions. The 
dominant local sources contributing to the ambient 
background are attributed to the emissions from 
several manufacturing factories (car industry; busi- 
ness parks; dockside) in the area and vehicular 
traffic emissions originating from two major free- 
ways surrounding the monitoring site from East to 
West (in a semicircle) and to a lesser degree lo- 
cal traffic on nearby roads. The freeways are lo- 
cated about l-2A;m away from the sampling site and 
carry, in general, high traffic volume throughout the 
day with the traffic peak hours between 07:00-10:00 
and 16:00-20:00. In terms of local biogenic sources; 
there is a park area about 1km SSW from the sam- 
pling site. The topography in the area is open and 
flat. 

3.3. Meteorological conditions 

Representative meteorological data have been ob- 
tained from three meteorological stations (Bureau 



of Meteorology sites) located 5-13A:to North, SW 
and SE from the sampling site. The air temperature 
and RH were in the ranges of 1-23°C & 36-100% 
with mean values of 11°C and 80%, respectively. 
Rainfall for the period was approximately 37mm 
per month, with an average of 12.5 rainy days per 
month. The prevalent wind conditions for the du- 
ration of the measuring campaign were WSW with 
average wind speeds of 2-4m/s. 

3.4- Data collection and processing 

Data were collected over 28 days between 
27/04/2009 and 23/06/2009. Approximately 500 
hours and 340 hours of data (5 seconds readings) 
were measured by FLAPS-III and APS, respec- 
tively. About 140 hours of data were measured 
simultaneously by both instruments. The data 
obtained were screened for outliers and anoma- 
lies, processed and stored in a custom built SQL 
database. A set of programming tools allowing data 
manipulation, retrieval and averaging was devel- 
oped using MATLAB. Exploratory and statistical 
analysis was performed using MATLAB and S-Plus 
software packages. 



4. Data Inspection 

Before discussing our turbulent mixing and dis- 
persion model, it is important to investigate some 
of the properties of the experimental measurements 
acquired throughout the campaign and validate the 
data set against similar studies. Simple statistics 
for APS and FLAPS-III data are presented in Ta- 
bles [T] and [21 respectively. 

4.1. APS Data 

The median values of particle number concentra- 
tion (dN) and mass concentration (dM) in the 1- 
10fj,m size range measured by APS were 1.5^/cm^ 
and 0.0067713/771"^, respectively. Median values are 
presented throughout this paper as a better esti- 
mate (compared to mean) due to skewed distr ibu- 
tion o f measured data d i stribu tion (Ho and Soencel 
()l998l ): [Morawska et all ^199^ ). 

A relatively low number concentration of par- 
ticles is expected, as the majority of particles in 
an urban/industrial type of environment is dom- 
inated by combustion sources (e.g. traffic emis- 
sions); generat ing predominan t ly fine , submicrom- 
eter particles (jJamriska et al ] (l2008h ; iMeiia et al 
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Parameter Measured 


Mean 


STD 


Median 


IQR 


Cone dN [1-lO^m] {#/cm^) 


2.189 


1.970 


1.500 


0.840 


CMD dN [l-lO^im] (pm) 


1.500 


0.101 


1.486 


0.107 


Cone dM [l-lO/im] {mg/m?) 


0.009 


0.007 


0.006 


0.004 


CMD dM [l-lO^im] ipm) 


3.122 


1.115 


2.839 


0.533 



Table 1: Statistics for aerosol data measured by APS (n ~ 245374; total sampling time ~ 340 hours). Particle mass concentra- 
tion dM was estimated from volume concentration d V and a particle density o{ ^ Ig jew? . dN and dM denote particle number 
and mass concentration; CMD denotes particles count median diameter. 



(|2fl07l) 1 which are outside of the size range consid- 
ered in this study. 

Particle mass concentration in the l-lO/im size 
range (dM[l-10/im]) derived from APS number 
data is also relatively low, reflecting the dominance 
of fine particles. An indirect comparison of the 
measured dM result s with literature data indicates 
a good agreement. iKevwood et al charac- 
terised PM concentrations in six Australian cities, 
reporting a mean value of PMi/PMio ratio of 0.45; 
i.e. approximately 55% of PMiq mass concentra- 
tion is attributed to particles in the 1-10/im size 
range. Using this estimat e and an ambient par - 



ticle densi t y of l .Gq/cTO (iMorawska et al.l (jl99 



Pitz et al.l (|2003l) ) the estimates for the PMio de- 
rived from our data are 0.026m(7/m'^ (mean) and 
0.017mf;/m'^ (median). These results are in good 
agreement with the annual mean PMio value for 
Melbourne of 0.023 ib 0.009mg/m^ reported by 
Kevwood et al. I (ll999l) as well as EPA data reported 
for the monitoring time period (0. 020-0. 02 4m(7 / rv?^ 
( Environmental Protection Agencv (|2009l )). 

In terms of particle size distribution, the corre- 
sponding median diameter for particle number and 
m,ass distribution was 1.5/irn and 2.8fim, respec- 
tively. The results presented are based on the anal- 
ysis of the available data for particles in the 1-10/im 
size range. While there is limited information or 
data on the size distribution statistics of this par- 
ticular size range, our results are in a relatively good 
agreem ent with findings reported by [Huffman et aT 
(I2OIOI) . 



In order to discern emission sources dominating 
the ambient background during different time pe- 
riods, the results were grouped into three time in- 
tervals: day (06:00-20:00), night (20:00-06:00) and 
24h (00:00-24:00) measurements. The statistics for 
these groups in the form of boxplots is presented in 
Figure [T] 

Visual observation of the data presented in- 
dicates that the mean and median values for 
the day, night and 24 hours data sets are com- 



a b 
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Figure 1: Box plot of particle (a) number {#/cm'^) and 
(b) mass (mg/m^) concentration measured by APS in 
the size range 1-10/im during day (06:00-20:00, nnay ~ 
153906); night (20:00-06:00), njv^ght ~ 91468); and 24 hour, 
(ii24(iotirs ^ 245374) periods. Mean value displayed as (□); 
median (centre line); 25% and 75% bottom and top of rectan- 
gle respectively; whiskers set to ±1.93(7 for display purposes; 
outliers (+). 



parable for both particle number and particle 
mass concentration, however statistical analysis 
of the medians difference (Mann- Whitney two 
sample rank-sum test) reveals a significant dif- 
ference between the medians of all three groups 
(at the 95% confidence level). There seems to 
be less fluctuation (i.e. spread of the data as 
indicated by IQR) in day data as compared to 
night data. This is likely associated with traffic 
emissions dominating the a mbient backgrouii d 
during day time (06:00-20 :00) (Mena et al. ('2007|); 



Bureau of Transport and R egional Economics [BTRE] 
(l2007l) 'l. 

To evaluate further the daily variation in APS 
number and mass concentrations, all available data 
were grouped and averaged over 30 minute time 
intervals. The averaged time series presented in 
Figure [5] show a complex multi peak profile with 
two broad dominant modes (6:00-13:00) and (16:00- 
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21:00) likely to be associated with an increased con- 
tribution from traffic during the traffic rush hours. 
In the morning, maximum number concentration 
was observed at about 09:00 2^ /cm?) while in 
the afternoon/evening mode the maximum was ob- 
served at about 20:00 (~ 1.7#/cm'^). There seems 
to be also relatively high concentration at night be- 
tween 23:00-03:00, which could be associated with 
the emissions from local industrial sources operat- 
ing 24/7 (e.g. dockside activities). 




4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 

Hour 



Figure 2; 30 minute average (Median it /QiJ/v^) of num- 
ber (dN, solid line) (#/cm^) and mass (dM, dashed line) 
{m,g/m^) concentrations measured by APS in the size range 
1-10/jm. (Number of data points n = 245374; 5 second read- 
ings). 

To investigate the trends further (in an attempt 
to unmask the effect of local emission sources) , the 
daily averages of particle number and mass concen- 
tration were normalised by the night average (i.e. 
median night values of dN and dM, respectively). 

Perhaps not surprisingly, the normalised concen- 
trations show similar trends as observed in Figured] 
with two dominant modes observed in the morning 
and late afternoon/evening. The concentration lev- 
els reached up to 1.4 times the night level between 
09:00-11:00 and up to 1.2 times at around 20:00 
(dN). Relatively low concentrations were observed 
for 03:00-05:00 and 13:00-15:00 time periods. The 
fluctuations could be associated with changes in 
emission sources (mainly traffic density) contribut- 
ing to the measured background and changes in at- 
mospheric conditions driving the aerosol dispersion 
(wind, convection, etc.). 

A similar trend as presented in Figure [5] was 
reported in a stu dy conducted in Brisbane by 
Meiia et al. ( 200?! ) for a comparable type of sam- 



pling environment (emission sources, location and 
topography). The authors showed a strong depen- 
dency of submicrometer particle concentration on 
traffic emissions (traffic density; type of cars) and 
wind direction. A review of the dependencies of am- 
bient background on the m eteorological conditio ns 
is presented for example in [Harrison et al.l (j2004l) . 



4.2. FLAPS-III Data 

The statistics for FLAPS-III data measured over 
a period of 500 hours is presented in Table [5J 
The SS counts represent the total aerosol parti- 
cles (TAPs) and FL counts the fluorescent biolog- 
ical particles (FBAPs) in the 1-lO^m size range 
measured over 5 second sampling intervals. FL1[2- 
31] and FL2[2-31] correspond to integrated particle 
count in bins [2-31] for the gated fluorescent light 
intensity as discussed previously. 

The median values for FLl/SS and FL2/SS ra- 
tio were 3.9% (FLl) and 5.1% (FL2), respectively. 
The mean values were 5.3% (FLl) and 7.2% (FL2). 
Correlation analysis of FLl and FL2 measurements 
(5 second readings) showed high correlation (R'^ ^ 
99%) indicating that fluorescence (from both FLl 
and FL2 channels) mostly originated from the same 
particle. Based on that assumption our results in- 
dicate that in terms of particle number, the FBAPs 
represent about 4-5% (median) of TAPs. 

Correlation between the fluorescent (FLl; FL2) 
and total (SS) count rate was also high (0.82 < 
< 0.84), which indicates a relatively high asso- 
ciation between the FL and SS data. 

Fluorescence in channel FL2 (emission wave- 
length ~ 580nm) was stronger compared to FLl 
(emission wavelength ^ 430nrn) with the FL1/FL2 
count ratio 0.79 (mean) and 0.75 (median). The 
effect could be associated with differences in fluo- 
rescence originating from different molecules (e.g. 
NADH and riboflavine) and diff e rent organisms 



( Agranovski and Ristovskil (j200 

2onnl) i. 



Brosseau et al. 



The presented results for total (SS count) and 
FBAPs are relative values; dependent on the phys- 
ical properties of sampled material, size distribu- 
tion of aerosols and FLAPS-III sampling character- 
istics. Using the absolute quantification of particle 
concentration capability of APS (TAPs), and the 
method described in Section [^31 concentration es- 
timates from FLAPS-III counts can be calculated. 
Although the instruments use different measuring 
techniques (APS: Time-of-flight, FLAPS-III: light 
scattering), a fairly acceptable correlation value of 
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Parameter Measured 


Mean 


STD 


Median 


IQR 


SS Total Count 


4376 


4754 


2849 


2684 


FL1[2-31]/SS Count Ratio 


0.053 


0.054 


0.039 


0.031 


FL2[2-31]/SS Count Ratio 


0.072 


0.071 


0.051 


0.044 


FL1[2-31]/FL2[2-31] Count Ratio 


0.788 


0.537 


0.745 


0.087 


FLl Total Concentration {^/cm?) 


0.046 


0.076 


0.023 


0.025 


FL2 Total Concentration [^/cw?) 


0.062 


0.108 


0.030 


0.033 



Table 2: Statistics for aerosol data measured by FLAPS-III. Analysis of SS Total Count and FL/SS Count Ratios include all 
data {N = 359571; total sampling time ~ 500 hours) measured in 5 seconds; Results for FLl and FL2 Total concentration 
were derived from FLAPS-III and APS paired data(Af = 99999). 



69.3% between the instruments was achieved over 
the measuring period. 

The measured fluorescent particle concentrations 
(mean values of 0.046 (FLl) and 0.062 (FL2) 
if/cm^) are in good agreement with the literature. 
[Huffman et"all ((2010) for comparison, reported the 
mean number concentration of FBAPs in the 1- 
10/ZTO size range of 0.03 (#/cm"^) - measured in a 
similar type of ambient environment. 

Further exploration of the fluorescent concentra- 
tion dependencies and behaviour was not the focus 
of this work and will be presented in a follow up 
publication. Since the measured FLAPS-III data 
will be used in relative terms as a ratio of FL/SS 
or normalised SS count, the terms count and con- 
centration will be used interchangeably throughout 
the following sections. 

Similar to the APS data, the FL/SS ratios de- 
rived from FLAPS-III data were analysed for dif- 
ferent time intervals (day, night and 24 hours) 
as presented in Figure [3l Several observations 
can be made from this data: the FL2/SS values 
were consistently higher compared to FLl/SS val- 
ues for both (day and night) datasets; the fluores- 
cent content was higher during the day compared 

to night ((FL/SS)Day > (FL/SS)Night); (FL/SS)Day 

values showed a narrower spread (IQR) compared 
to night results and were shifted towards the lower 
values (smaller fraction of fluorescence). Analy- 
sis of the medians difference (Mann- Whitney two 
sample rank-sum test) showed statistically signifi- 
cant difference between the medians of FLl/SS and 
FL2/SS for day, night and 24h time periods (at the 
95% confidence level). 

The daily variation in fluorescent particle load- 
ing calculated as 30 minute averages is presented 
in Figure |31 The fraction of fluorescent particles 
varied throughout the day between 3-7% (FLl/SS) 
and 4-9% (FL2/SS). The variation has diurnal char- 
acter with the most dominant peaks appearing at 




Day Night 24 Hours 



Figure 3: Box plot of (a) FLl/SS and (b) FL2/SS count ratio 
measured by FLAPS-III in the size range 1-lO^m during day 
(06:00-20:00; tioay ~ 216188), night (20:00-06:00; ^Night ~ 
143383) and 24 hour {N24,hours ~ 359571) periods. Mean 
value displayed as (□); median (centre line); 25% and 75% 
bottom and top of rectangle respectively; whiskers set to 
±1.58(7 for display purposes; outliers (-I-). 



around 11:00 and 15:00. The fluorescent particle 
content showed a steady increase starting at about 
5:00 until 9:00, with a sharp increase peaking at 
10:30. A second, less pronounced maximum can be 
observed at about 15:00. 

To discern the dependencies and behavior of the 
24h daily data, the daily variation of FL and SS 
data were normalised by the night values as pre- 
sented in Figure [S] It can be seen that the SS data 
shows several peaks throughout the day (the most 
dominant being between 10:00-12:00) while FLl 
and FL2 signals exhibit two broad distinct modes 
(peaks). Similar to the SS data, the most domi- 
nant peaks for FLl and FL2 are observed at about 
10:00. The FL levels start to increase at about 
05:00 which could be associated with the emission 
from biogenic sources linked to sunrise with a fur- 



7 




12 3 4 5 6 7 8 9 101112131415161718192021222324 
Hour 



Figure 4: FLl/SS (solid line) and FL2/SS (dashed line) Ra- 
tio's, 30 minute median averages, IQR/^/n shown as error 
bars. 



ther steady increase between 08:00-09:00 and then 
a very rapid climax with a peak at 10:30. We spec- 
ulate that the dominant mode during 10:00-12:00 
period could be associated with the presence of bi- 
ological material originating from different biogenic 
sources (e.g. marine environment) emerging as a re- 
sult of changed environmental conditions (e.g. wind 
direction). Other factors influencing FL/SS fluc- 
tuations could be attributed to different biological 
material c omposition with d ifferent strengths of flu- 



orescence (Cozetal 



( 2010f) '). Investigation of the 
composition of the biological materials and their 
source apportionment was beyond the scope of this 
study. 
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Figure 5: 30 minute averages (Median) for FLAPS-III parti- 
cle count normalised by the median night count for total (SS 
[solid line]) and Fluorescent (FLl [dashed line]; FL2 [dotted 
line]) particles in the 1-10/^m size range. Median night val- 
ues of 2430 for SS, 72 for FLl and 91 for FL2. 359570 data 
points. 5 second samples. 



An indirect validation of our results was obtained 
through comparison with the results of other stud- 
ies focused on biobackground characterisation using 
o ptical and bio l ogical methods as outlined below. 



Boulet et al.l (|1996l ). based on a review of previ- 



ous studies in different locations, showed that under 
normal, steady state background atmospheric con- 
ditions the fraction of fluorescent particles rarely 
exceeded 5%. 

IVIonitoring of fluorescence in ambient air mea- 
sured in Calgary CA, where relatively high lev- 
els of human and vegetative contributions could 
be expected, the fluorescent fraction of ambient 
aerosol in the 2.5-lO^m range rarely exceeded 5% 
regar dless of the tot a l part icle number concentra- 
ti on dHo and Speiicel (|l998[) ). 

Coz et al. ( 2010t ) reported that PB APs measured 



in an urban area in the northeastern Unites States 
contributed to PM2.5 on average 6.9 ± 5.4% during 
summer and 3.3 ± 1.4% during winter. The high 
variability during summer could be attributed to 
short term variations in both the load and compo- 
sition of bacterial aerosols. The small variability 
found during winter suggest the existence of a rel- 
atively constant background level of PM2.5 PBAPs 
mass of 2-5%. 

In a study v e ry sim ilar to the presented work. 



Huffman et al 



(l2010l) characterised the bioback- 
ground in Mainz, Germany in a semi- urban envi- 
ronment. Using similar experimental and sampling 
setups (FLAPS-III sampling from a building roof 
close to the CBD) over several months, the authors 
reported about 4% of FBAPs in terms of particle 
number and about 20% in terms of particle mass. 

In summary, the measured data are in a good 
agreement with results reported in literature in 
terms of TAP and FBAPs concentration, fluores- 
cent content (FL/SS) and their fluctuation during 
the day. These results imply validity of the col- 
lected data set and its suitability for further explo- 
ration using advanced statistical methods aiming to 
develop superior alarm trigger algorithms. 



5. Model of Biobackground 

Particles in the biobackground are in continuous 
motion (with random displacements) and are be- 
ing advected by turbulent flow in the atmospheric 
surface layer. Due to their size (l-lO/zm) the in- 
verse effect of these particles on flow is negligible 
and their kinematics can be described by the mod- 
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els of 'passive scalar ' turbulence (or turbul ent La- 
grangian transport) iFalkovich et al. ( 2001 ). This 
means that one can introduce the instantaneous 
concentration of the particles 9{r,t) that obeys a 
stochastic equation of concentration transport (the 
stochasticity of the transport equation is a result of 
the randomness of the underlying turbulent velocity 
field). In other words, the theory of scalar turbu- 
lence provides a consistent framework for deriva- 
tion of the statistical properties of concentration 
0{r, t) (i.e. moments, PDFs and correlations) for 
given statistics of velocity. Remarkably, some im- 
portant properties of particle statistics experimen- 
tally observed in complex environmental settings 
(boundary layer flow, strong anisotropy, buoyancy 
etc.) emerge from the simplest (but still justifi- 
able) models of tur bulent flow i.e. white- no ise or the 
Krachnan model (jPalkovich et all (|200l[ )). From 
this perspective we anticipated that certain statis- 
tical characteristics of the biobackground should 
also reflect the fundamental properties of under- 
lying mixing processes, irrespective of the nature 
(biological/non-biological) of the aerosol particles. 
This was the main rationale behind our approach 
to the exploration of our experimental data and the 
development of a statistical framework as outlined 
below. 

One of the most interesting predictions of the 
theory of scalar turbulence applied to transport in 
the atmospheric surface layer is the conclusion that 
all statistical moments of the scalar concentration 
{6"{r,t)) are expected to follow the same proflle, 
regardless of the value of n. This allows one to 
write the following simple scaling la ws for the con- 
centra ti on moments (for detafls see Yee an d Chan 
( 1997t). iLebedev and Turitsvnl (|2004l ). lYeei (,20091 
Skvortsov and Yed (|201llv r 



n> 2 



and 



(1) 



(2) 



The pre-factor ^„ in this expression has been 
theoretically evaluated for two extreme cases. For 
the case of a perfectly mixed tracer (which usu- 
ally corresponds to locations far away from a tracer 
source) 1^,7 — 1; while for the opposite, highly 
intermittent case (near the source region) ^+ = 
r(rt -I- l)/2"~^; whe r e r(-) is the Gamma function 
( Skvortsov and 'Ye^ (|201lh . Hence, m real situa- 



tions we expect that particle concentration statis- 
tics should follow the scaling law ([Ij with the 
pre-factor ^„ located between these two estimates 

(C < e« < a)- 

We found these predictions are in good agree- 
ment with our extensive observational data (see 
Figures [H and [7|) . The statistical moments of par- 
ticle concentrations {9'^{r,t)) were calculated from 
time series recorded by the FLAPS-III (SS, FLl, 
FL2) and APS instruments at the same location, so 
d(r, t) = 9{t) in our case. The scaling laws ([T]) corre- 
sponding to the experimental data from FLAPS-III 
for n — {3,4} are depicted in Figure [5] and com- 
pared with the theoretical predictions ([!]) and ([5]). 

The non-biological contribution (top plot in Fig- 
ure [5]) comes from the SS channel of FLAPS-III 
and the biological (bottom plot in Figure (5]) from 
the summation of the FLl and FL2 channels. Since 
FLl and FL2 are strongly correlated the aggregated 
FL1+FL2 metric captures all relevant data and re- 
moves the need to analyse both channels separately. 
APS concentrations follow the same trends and are 
therefore not shown in the following analysis. 

The experimental data follows the scaling laws 
(HI) reasonably well within a wide range of measured 
concentrations (see Figure [3] and Table [2|) with the 
pre-factor value ^„ extensively occurring in the ex- 
pected range ^~ < < ■ For the n = 2 case 
^ , a similar plot is presented in Figure [T) Since 
pre-factor ^2 is not deflned in ([2), a value of 3.5 was 
assigned for better visual appearance on the theo- 
retical fit. From Figure [HI we can see that when ap- 
plied to the biobackground the scaling laws ([1]) , ^ 
become less accurate for higher order moments (i.e. 
higher n) from the presence of outliers above the S^f 
bound on the bottom plot of Figure [6l This is sim- 
ilar to t he case of tracer di s persio n from localised 
sources (|Yee and Skvortsovl (|201ll )'l. Nevertheless, 
with the two theoretical values of pre-factor they 
always provide a valuable tool for the estimation of 
the statistical moments of the aerosol particle con- 
centration in the atmospheric surface layer. 

To validate our scaling conjecture further and to 
prove its independency of the particles' origin we 
present the scaling law ^ of biological and non- 
biological sources in Figure [T) We observe that the 
scaling law ([2]) of the statistical moments of the 
biobackground concentration are in excellent agree- 
ment with the prediction of the scalar turbulence 
theory, i.e. they are mostly determined by the kine- 
matics of the turbulent surface layer. 

It is worth stressing that the scaling laws ([T]) and 
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Figure 6: Scaling laws for statistical moments of particle 
concentration. Top plot represents non-biological material 
(SS) and the bottom biological (FL1-I-FL2). Bottom and 
top data series on each plot correspond to moments n = 3 
(o) and n = 4 (□) respectively. Lines correspond to the 
theoretical relationship 1^ with two different values of pre- 
factor ^„ = 5^ (dashed line), 5„ = (solid line), see text 
for details. 



([21) manifest a general property of turbulent mixing 
and are not a consequence of any particular func- 
tional form of concentration PDF p{9). Contrarily, 
any distinct PDF form should comply with these 
scaling laws and this compliance is one of the im- 
portant criterion for the evaluation of appropriate 
candidate model for p{9). 

Amongst the different functional PDF forms 
of tracer concentration evaluated in studies of 
turbulent tracer dispersion, the Gamma form has 
recently emerged as one that adequately captures 
the physics of turbulent transport and provides 
a good match with a variety of experimental 



Figure 7: Similarity of scaling laws for statistical moments 
of particle concentration for biological (FLl: □, FL2: o) and 
non-biological (SS: 0) particles, see text for details. The 
dashed line corresponds to the theoretical prediction JJJl, 
with a scaling constant of 3.5 for visual appearance only. 



data ( Yee and ChanI (1997); Skvortsov and Ye3 
(|201l[ )). This model (see the Gamma form of 
(j31)) can be analytically derived using a phe- 
nomenological analogy between the mixing and 
convolution processes by explicit modeling of 
the deformation of trac er blobs by random 
stretc hing a nd folding ([Duplat and Villermaux 
( 20081) : IVillermaux and Duplatl (|2003l ): 

Venaille and Sommeria ( 20081 )). Indeed, the 
Gamma PDF complies with the scaling laws 
of the concentration moments ((Il),(l2]) with 
the pre-factor ^„ residing within the range 
£^ < £ra < d e scribe d above (for details see 
Yee and Skvortsov! ( 2011 )). Whilst the Gamma 
distribution is defined with two parameters 
(scale 6 and shape fc), it can al so be written 

in a "single-pa r a.meter" form dYee aiid ChanI 

(I1997I): 1^1 2009'): 'Duplat an d Villermauxl (l2008h : 
Villermaux a nd Duplat (2003)): 



p(x^0/(0)) =^x'^-^exp(-A;x), 



(3) 



where k is the only fitting parameter of the model 
and can be related to the statistical moments by 



(|Yee and Skvortsov! (|201lf )) 



1 
k 



{0? 



- 1. 



(4) 



From these two equations we can see that our 
model of the biobackground is uniquely determined 
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by the first two moments (i.e. (9) and {0'^)), which is 
a well known property of the Gamma distribution. 





Figure 8: CDFs of non-biological (SS, top plot) and biologi- 
cal (FL1-I-FL2, bottom plot) particle concentrations (o) with 
intermittency correction ((SJl. Gamma (solid line) and Expo- 
nential (dashed line) models are also shown. Inserts show 
the Exceedance Distribution Function (1-CDF) of the same 
data and models, see text for details. 

To verify the Gamma PDF model agrees with our 
experimental data, we compare it with the statisti- 
cal properties of our observations (i.e. time series of 
particle concentrations from FLAPS-III). In order 
to do this; (9) is calculated for the bio- and non 
bio- particles, then histograms of our time series of 
concentration in terms oi x = 9 / {9) were produced. 
These histograms (or empirical Cumulative Distri- 
bution Functions (CDFs)) are compared with the 
CDF corresponding to the functional form of 
In such an approach, k was the only parameter to 
fit from observations. 

The results of this comparison are depicted in 
Figure [8] (non-biological particles top; biological 



bottom). The tight fit shows that this model can 
be reliably used for the parametrisation of the sta- 
tistical properties of ambient air (both biological 
and non-biological contributions). The estimations 
of parameter k were quite stable leading to values 

= 1.1 for non-biological and A: = 1.8 for biological 
particles. For non-biological particles the (mathe- 
matically) simpler Exponential PDF (i.e. ^ with 
fc = 1) can be used as an approximate model for 
our data. This distribution has a k value of 1.1 and 
therefore a relatively small deviation of 0.1 from the 
k = 1 required for the Exponential form. The sig- 
nificant advantage of the Exponential model (even 
taking into account its inferior data fit compared to 
Gamma, see both plots on Figure [8]) is that it re- 
quires fewer parameters to estimate and to fit. For 
instance; the Exponential PDF model can be used 
for rapid calibration of a detection system when 
only a limited set of data for the monitored envi- 
ronment has been collected. The biological contri- 
bution (bottom plot of Figure [5]) is described quite 
poorly by the Exponential PDF. This contribution 
(the presence of biological material) can therefore 
be easily identified just by monitoring the value of 
k. Investigation of this dependency will be the first 
parameter of interest in our future publication on 
biodetection algorithms. 

The Gamma model ^ does not take into ac- 
count the intermittency of a particle distribution 
(i.e. a finite probability that the particle concentra- 
tion is exactly zero for a certain period of time). 
This model is therefore only valid for perfectly 
mixing tracer flows - which rarely occur natu- 
rally. In order to estimate the effect of inter- 
mittency we appropriately modifled our concentra- 
tion PDFs by adding an inter i nitten c y co r rectio n 
(for details see lYee and ChanI (Il997l) : PYe^ (flooi); 
Skvortsov and Yed (12011^ " and references therein) 

(5) 



p{9) = {l-^)5{9)+-ip{9), 



where 5{-) is the delta function, 7 is the inter- 
mittency factor (a degree of mixing with 7=1 
corresponding to the perfect mixing case), p(0) is 
the "regular" PDF (in our case ([3])). The estima- 
tion of parameter 7 in the model ([5]) from concen- 
tration time series is well e stablished in the stud - 
ies of turb ulent disp ersion ( Yee and C hanI ( 1997 ): 
Ye^ (|2009t ): .Skvort sov and Yee (2011.) and refer- 
ences therein), so we omit all technical details and 
present only the main results. We found that the 
effect of 7 is not signiflcant in our data, i.e. 7 was 
always close to unity. A stable estimation from 
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both FLAPS-III and APS yield a consistent value 
7 ~ 0.97, which is exemplified in the small positive 
shift in the p(x) -intercept from in Figure |S1 Nev- 
ertheless, since the value of 7 is also determined by 
the proximity to a particle's source, we expect that 
in some cases the intermittency co rrection in the 
model ( 151 inay be come appreciable (|Yee and Chan 
(|l997h :]^ (l2009h 'l. 



If this is the case, 7 will be the second parameter 
of interest in our biodetection algorithm. This as- 
sumption also leads to two important conclusions. 
Firstly, there should be a significant similarity be- 
tween bio- and non bio-aerosol characteristics mea- 
sured at the same location (since they undergo the 
same mixing process). Secondly, these characteris- 
tics should clearly comply with the physics of tur- 
bulent mixing. As illustrative examples we mention 
here a universal relationship between statistical mo- 
ments of particle concentration (a so-called scaling 
law) and the explicit functional form for the prob- 
ability distribution function of particle concentra- 
tion (the Gamma function). It is worth remark- 
ing in this context, that the scaling law for statis- 
tical moments is more general and is not bound to 
a particular functional form o f concentration PDF 
(iLebedev and Turitsvnl t004)). 



6. Conclusions and Future Work 

In summary, we quantified the content of am- 
bient background for an urban/industrial type of 
environment close to the Melbourne CBD. Concen- 
tration of total and biological aerosols in the 1- 
10/im size range was monitored over a period of 
two months and the results of exploratory analyses 
for day, night and daily (24 hour) fiuctuations are 
presented. 

The measures of total particle number and mass 
concentration are in good agreement with previous 
results obtained for the Melbourne area. Daily pro- 
files (30 minute averages) shows complex behaviour 
with two dominant modes observed in the morn- 
ing and afternoon/evening corresponding to traffic 
peaks hours. The results indicate that the ambi- 
ent background loading is dominated by vehicular 
traffic emissions. 

The fraction of fluorescent particles in total was 
in the range of 3-9% (30 minute average) which is in 
good agreement with other studies reported in the 
literature. Concentration levels of biological mate- 
rial showed diurnal variation which could be asso- 
ciated with an increase in emissions and changes of 



biogenic sources. The levels increased up to four 
times during the midday period as compared to the 
night mean values. O ur results are in good agree- 
ment with the study [Huffman et al. con- 
ducted in a similar type of ambient environment. 

A universal model of the biobackground has been 
produced from this data set that provides a rigorous 
framework for the development of novel algorithms 
for bio-aerosol detection/characterisation and can 
help to both improve and optimise some existing 
algorithms used in operational systems. 

These refinements may be implemented by ini- 
tially using the scaling laws (IT|) and ([2]) , which pro- 
vide an efficient way to estimate any statistical mo- 
ment from a limited data set. For instance, mo- 
ments can be obtained even if only a few concentra- 
tion levels are known (measured) or measurements 
are of poor quality. This implies that by employing 
relationship ([TJ we can continuously infer the mix- 
ing state of the biobackground, i.e. how close its dis- 
tribution is to the equilibrium state. In effect, any 
deviation from the equilibrium scaling given by law 
(IT|) would indicate the proximity of a bio-aerosol re- 
lease. Secondly, the specific PDF form the bioback- 
ground (i.e. Gamma PDF) enables a rigorous as- 
signment of a well-defined measure (based on like- 
lihood or risk) that is associated with any detected 
spike in bio- aerosol measurement. An informative 
decision regarding the relevance of this spike can 
then be made according to the current operational 
context (e.g. it may provide a consistent method- 
ology to reduce the false alarm rate), resulting in 
rigorous risk mitigating strategies. Lastly, the deep 
analogy between bio- and non-bio particle distribu- 
tions emerging from the physics of turbulent mix- 
ing allows significant simplification of any calibra- 
tion process for deployable systems (including the 
capability for rapid self-calibration), by extracting 
additional statistical information from monitoring 
the non-biological aerosol background. 

A monitoring system such as this would en- 
able conscious operational decisions to be employed 
which can help to develop more robust instrumen- 
tation algorithms for the detection of biothreats in 
operational settings, rigorously evaluate the risk as- 
sociated with miss and false detection as well as to 
propose relevant mitigation strategies in response 
to the detected bio-release. Specifics concerning 
this framework will be elaborated upon in future 
publications. 
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